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Abstract. We study the two-pion exchange potential at next-to-next-to-leading order in chiral effective 
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Z > 2 as compared to the calculation utilizing dimensional regularization. 
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' 1 Introduction 

' Effective field theory (EFT) has become a standard tool 
, for analyzing the chiral structure of Quantum Chromody- 
' namics (QCD) at low energies, where the perturbative ex- 
[ pansion in powers of the coupling constant cannot be used. 
■ The chiral symmetry of QCD is spontaneously broken and 
' the corresponding Goldstone bosons can be identified with 
the pions, if one considers the two flavor sector of the up 
' and down quarks as done here. It is a general property 
[ of Goldstone bosons that their interactions become weak 
• for small momenta. Chiral Perturbation Theory (CHPT) 
' is the effective field theory of the Standard Model which 
allows to describe the interactions of pions and between pi- 
ons and matter fields (nucleons, p-mesons, Zi-resonances, 
. . .) in a systematic way. This is achieved via an expansion 
of scattering amplitudes and transition currents in powers 
of a low-momentum scale Q associated with small exter- 
nal momenta and with the pion (light quark) mass. Here, 
small means with respect to the scale of chiral symmetry 
breaking of the order of 1 GeV. Pion loops are naturally 
incorporated and all corresponding ultraviolet divergences 
can be absorbed at each fixed order in the chiral expan- 
sion by counter terms of the most general chiral invariant 
Lagrangian. 

This perturbative scheme works well in the pion and 
pion-nucleon sector, where the interaction vanishes for 
vanishing external momenta in the chiral limit, for some 
early reviews see e.g. 1,2 ,3 ,4 . The situation in the few- 
nucleon sector is much more complicated. The main dif- 
ficulty in the direct application of the standard methods 
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of CHPT to the two-nucleon (2N) system is due to the 
non-perturbative aspects of the problem, the unnaturally 
large S-wave scattering lengths and the shallow nuclear 
bound states. One possible way to deal with this difficulty 
has been suggested by Weinberg, who proposed to apply 
CHPT to the kernel of the corresponding integral equa- 
tion for the scattering amplitude, which can be viewed as 
an effective nucleon-nucleon (NN) potential [321 ■ 

The first quantitative realization of the above idea has 
been carried out by Ordofiez and co-workers, who de- 
rived an (energy-dependent) NN potential and performed 
a numerical analysis of the 2N system [7j. The energy- 
independent representation of the chiral NN potential, 
which can be applied much easier in few-nucleon calcu- 
lations, has been derived in "WW^IOK At leading order 
(LO) in the chiral expansion the potential is given by the 
well established one-pion exchange (OPE) and two con- 
tact forces without derivatives. At next-to-leading order 
(NLO) OPE gets renormalized and the leading two-pion 
exchange (TPE) diagrams as well as seven more contact 
operators appear. At NNLO, one has to include sublead- 
ing TPE with one insertion of dimension two pion-nucleon 
vertices (the corresponding low-energy constants (LECs) 
are denoted by ci_3,4, we adhere to the notation of Ref. Pj). 
Notice that no new contact forces contribute at this order. 
While the pion exchanges are governed by the underlying 
chiral symmetry, the contact forces represent our igno- 
rance of the short-range physics and are not restricted by 
chiral symmetry. 

The LECs 01^34 enter the expressions for the pion- 
nucleon (ttN) scattering amplitude at subleading (Q^) and 
higher orders and thus can be determined from the ttN 
data. From the analysis ^Jj one gets: ci = —0.64 , C3 = 
—3.90 , C4 = 2.25 (here and in what follows the values of 
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the Ci's are given in GeV ^). The values obtained from 
various analyses jlllll2lll8lll4lll5j are in the ranges: 
ci = -0.81 . . .-1.53 , C3 = -4.70 . . .-6.19 , ca = 3.25 . . . 4.12 
These bands are also consistent with expectations from 
resonance saturation [^. Notice that the numerical val- 
ues of C3 and C4 at both orders and are quite 
large, which can be partially explained by the fact that 
the LECs C3^4 are to a large extent saturated by the A~ 
excitation. This implies that a new and smaller scale, 
namely rriA — m ^ 293 MeV, enters the values of these 
constants in EFT without explicit Z\, see ^^l- At fourth 
order, these LECs get modified by quark mass dependent 
contributions and also are affected when electromag- 
netic corrections are included |2| ■ These modifications go 
beyond the accuracy of the calculations performed here 
but should be kept in mind if one wants to go to higher 
orders and/or systematically includes isospin violation. 

The large numerical values of the c^'s lead to dramatic 
consequences in few-nucleon systems |18) . The resulting 
subleading TPE correction calculated using dimensional 
regularization |^ (or equivalent schemes) turns out to be 
very strong already at intermediate distances r ~ 1 — 2 
fm. This could, in principle, be compensated by the cor- 
responding contact terms. However, such a compensation 
at NNLO is only possible in low partial waves (i.e. in S-, 
P-waves as well as for ei) since the contact terms do not 
contribute to D- and higher partial waves at this order. 
The D- and F-waves may therefore serve as a sensitive 
test of chiral TPE exchange, as suggested by Kaiser et al. 
in in] , since higher partial waves are strongly dominated by 
OPE and less sensitive to TPE (as it is known since long, 
for an early nonperturbative approach see jJH). The con- 
ventional scenario of nuclear forces represented by exist- 
ing OBE models and various phenomenological potentials 
suggests that the D- and higher partial wave NN inter- 
actions are weak enough to be treated perturbatively, see 
|18) . Clearly, in such a framework one can not describe 
the low partial waves that show strong nonperturbative 
effects. Under this assumption Kaiser et al. [H] applied 
chiral EFT to perform a parameter-free calculation for 
the neutron-proton {np) D- and higher partial waves and 
found rather poor convergence of the chiral expansion al- 
ready at surprisingly low energies. While the LO and NLO 
predictions, which correspond to OPE and to OPE accom- 
panied by the leading TPE corrections, already agree rea- 
sonably well with the Nijmegen phase shift analysis [20] 
(NPSA), taking into account subleading TPE at NNLO 
governed by the LECs 01,3,4 destroys that agreement and 
leads to strong deviations from the data for E\sh > 50 
MeV in D-waves and for Slab > 150 MeV in F-waves. 

In |2T] we have demonstrated that the NNLO poten- 
tial allows for a good description of the NN data, which is 
also visibly improved compared to the NLO results. Con- 
trary to Ref. [H], we did not perform a perturbative ex- 
pansion of the np T-matrix in high partial waves and cal- 
culated phase shifts by solving the Lippmann-Schwinger 
(LS) equation for the NN T-matrix. We found that taking 
the momentum space cut-off in the LS equation of the or- 
der of 1 GeV allows for a satisfactory description of all par- 



tial waves simultaneously. With such a large value of the 
cut-off, the isoscalar central TPE potential becomes al- 
ready so strongly attractive that unphysical deeply bound 
states appear in the D- and in the lower partial waves. 
Note that although such deeply bound states do not in- 
fluence NN observables at low energies, they might show 
up in another processes (like e.g. Nd 18 and wd [221 scat- 
tering). Since the potential is very strong (and attractive) 
and there are no counter terms according to the power 
counting, changing the value of the cut-off clearly leads 
to a strong variation of the ZJ-wave phase shifts. Higher 
order counter terms are needed in order to reduce the cut- 
off dependence of these observables and thus the problem 
with the slow convergence of the chiral expansion remains. 

Motivated by the known cancellation between the tttt 
and np exchanges, which has been observed in boson ex- 
change models of the nuclear force, we have constructed in 
Ref. [IHl the NNLO* version of the NN potential without 
spurious deeply bound states. To achieve that, we adopted 
values of the LECs 03,4, which are much smaller in mag- 
nitude than the ones obtained from ttN scattering and re- 
sult from subtracting the /i-contribution and fine tuning 
to NN observables: C3 = -1.15 GeV-\ a = 1.20 GeV-^. 
This also allowed for a fairly good description of the D- 
and higher partial waves. Accounting for subleading TPE 
leads to small corrections in most channels and the chiral 
expansion for NN scattering seems to converge. Certainly, 
the situation is still far from being satisfactory since the 
small values of LECs 03,4 are not compatible with ttN scat- 
tering. Notice that the large values of these LECs are also 
supported by recent determinations from pp and np par- 
tial wave analysis performed by the Nijmegen group j23l 

El- 

In the present work we explain the origin of the above 
mentioned problems and present a way to improve the 
convergence of the chiral expansion for the NN interac- 
tion. It allows to use the large values of c^'s consistent 
with ttN scattering. We argue that the unphysically strong 
attraction in the isoscalar central part of chiral TPE at 
NNLO resulting when calculated using dimensional (or 
equivalent) regularization is due to high-momentum com- 
ponents of exchanged pions, which cannot be properly 
treated in an EFT. Using a cut-off regularization instead 
of dimensional one and taking reasonable values for the 
momentum space cut-off allows to remove spurious short- 
distance physics associated with high-momentum inter- 
mediate states and to greatly improve the convergence 
of the chiral expansion. A similar idea has already been 
used a long time ago in the analysis of the octet baryon 
masses and the pion-nucleon sigma term [2S| and was re- 
cently applied to improve the convergence of the SU(3) 
baryon chiral perturbation theory ^26,27,28 . A critical 
discussion about the use of cut-off schemes is provided 
in Ref. [221 ■ Notice further that a finite momentum-space 
cut-off in chiral loops has been used to derive the (energy- 
dependent) expressions for TPE in Ref. [J. We also pro- 
pose a simple and convenient way to derive analytic ex- 
pressions for regularized TPE in the momentum space 
based on the spectral function representation. 
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Our manuscript is organized as follows. In section |21 
we describe our formalism and present the explicit ex- 
pressions for regularized TPE. In section |2| we apply the 
formalism to np D- and higher partial waves and compare 
the results with the ones obtained from dimensionally reg- 
ularized expressions. The summary and conclusions are 
given in sectional 



2 Formalism 

Energy-independent expressions for the chiral TPE at NLO 
and NNLO have been derived using different formalisms 
in Refs. 8 !) ID. The corresponding diagrams are shown 
symbolically in Fig. The last diagram in the first line 
(NLO) requires a special treatment in order to avoid dou- 
ble counting of the iterated OPE. Notice further that the 
first two graphs in the second line lead to vanishing con- 
tributions to the NN force. The explicit time-ordered ex- 
pressions for the potential can be found in jJUj. While 
OPE is of the order (Q/yl^)*', where refers to the chi- 
ral symmetry breaking scale of the order of the p-meson 
mass, TPE at NLO and NNLO provides corrections of or- 
ders {Q/Ay^Y s^iid {Q/A^)^, respectively. In the present 
work we will adopt the counting rule for the nucleon mass 
m ^ A^, which is required for Weinberg's power counting 
to be consistent (SjjT,. Therefore, we do not need to in- 
clude relativistic l/m-corrections (1/m^-corrections) to 
TPE (OPE) at the order considered. Notice that in the 
one-nucleon sector m is usually treated on a same foot- 
ing as A^. The same counting rule has also been used in 
Ref. inj. The TPE contributions at NLO are given by: 



dH (1^-q^) 



(27r)3 (ijj_|_ -I- u}_) 



1 , . f dH {uj+-uj^f 1 

4 ("^1 ■ '^2) 



8(2F,) 

2(2F^)'i y (27r)3 uj\ujt(uj^+ijj_) 
x((Ti.X2) (l^-q^ 



-6(a2-[qxl])(ai.[qxl]) 



(2.1) 



where Oi and 
the nucleon 



c.± = v/(^±iF 

reads: 



Xj are the spin- and isospin-matrices of 
q is the nucleon momentum transfer and 
The subleading TPE potential 



■4M2. 



V 



NNLO 



3gi 

16F4 7 (27r)3 
X f 8ci Ml - 



dH 12 -q^ 



x(ct2 



C3 (I'-q^ 
dH 



^'^^ ■ J (27r)3 LU+LU_ 

[qxl])x(ai.[qxl]). 



(2.2) 



We will show in the next section how these integrals can 
be evaluated using different regularization schemes. 



2.1 Dimensional versus cut-off regularization 

The integrals in Eqs. 1)2. l|l . (|2.2|) are ultraviolet divergent 
and thus need to be regularized. Applying dimensional 
regularization the TPE potential takes the form : 



VdR — VuR, non-pol. + ^DR^ pol. 



(2.3) 



where the non-polynomial and polynomial (in q ) parts at 
NLO F^LO^^^ p^^ and V^^% , respectively, read: 



T^NLO 
DR, non— pol. 



3847r2i^4 



+ q^{23g\-10g\-l) 



L{q)\AMl{5g\~4g\-l) 



48g\M^ 



L{q) (ai.q)(fT2-q; 



T/NLO 
•'DR, pol. 



647r2F4 

-((71 • 0-2) 



(2.4) 



+S: 



DR 



(cti • q) ((72 ■ q) - ((71 • (72) q" 



Here, we have set (7 = |q | and the logarithmic loop func- 
tion L{q) is given by 



q 2AU 



+ 4M2 . (2.5) 



Further 

= 



5^{(-23.UlO.Ul)ln^ + .4, 

(2.6) 



647r2Fi 



where fj, is the scale of dimensional regularization and the 
ai are polynomials in gA, whose precise form depends on 
the choice of subtraction. The non-polynomial part as well 
as all terms proportional to ln(A/^//.t), which are due to 
logarithmic divergences in Eq. (|2.1() , are unique and do not 
depend on the choice of subtraction. The /i-dependence 
of the V[^R°oi is compensated by the corresponding NLO 
counter terms of the form 



NLO 



+ Ciip) (c^i • q)(CT2 ■ q 



C^2(m) q^ (T1-T2) 

)+Ci{pL)q^ ((71 •CT2) 



(2.7) 



and the resulting renormalized NN potential does not de- 
pend on /i. Application of dimensional regularization to 
nuclear potentials is also discussed in |3()| . 
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NLO: 



y (1- 



1. (h 



NNLO: 



Fig. 1. Chiral TPE at NLO and NNLO. Heavy dots denote leading vertices from the chiral 
Lagrangian while solid rectangles correspond to subleading ones, which depend on the LECs £1,3,4. 



At NNLO one finds: 



T/NNLO 
DR, non— pol. 



3gi 



2M2(2ci - C3) - cag' 



^{2Ml + q^)A{q) 



9\ 



v. 



NNLO 



DR, pol. 



327rF4 

X (cti • q)(cr2 ■ q) - g^(cri • (T2 

( qDR ^ qDR «2\ I qDR / ^ \ 



(2. 



X 



((Ti • q) (cr2 ■ q) - (o-i • (T2) (t 



with the loop function A{q) 



. , , 1 q 
A\q) = — arctan — — — , 



(2.9) 



and 



5™ - - 



167rFi 



C3)M3 



(2.10) 



cDR _ gj 

3 - 327ri^i 



Notice that the integrals in Eq. H2.2|) are finite in dimen- 
sional regularization. 

Although dimensional regularization provides an easy 
and convenient regularization scheme, it is by no means 
the only possible one. One can equally well regularize the 
divergent integrals in Eqs. H2.HI . 12.2|l using a momentum 
space cut-off, i.e. by multiplying the corresponding inte- 
grands by the regulating function /a(0 = //i(|l|), with the 

properties /a(0 1; /a(0 0- This function /a(0 
should go to for large I quick enough in order that the 



regularized integrals exist. For cut-off regularized (CR) 
TPE one finds similar to Eq. ifT^ : 



VcR = VcR, non-pol. + VcR, pol. , 



(2.11) 



where 



Vr 



NLO 

CR, non — poL 



°p T/NLO 



DR., non— poL 



T/NNLO A^oo ,^NNLO 



CR, non— poL 



DR., non— poL ' 



(2.12) 



The polynomial pieces have the same structure as in Eqs 
(|2.8I) . where the quantities 5°^ (Sf^) should now be re 
placed by ^ (^f") given by 



5p^ = 



^-—l-l8M^{5g\ - 2g\) In ^ 

+a'iAf2 + /3iyl2| ^ (2.13) 



cCR _ 

^ 3847r2^^4 

cCR _ "^9 A 
^ 64^2^4 



In — + "3 



and 



CR 



cCR _ 

' ~ 16^Fi 



35^ 

{c3M^ + (3iA), 



(2.14) 



cCR _ 
D3 — 



9\ 



327ri^i 



(c4Af^ + ^5^) • 



Here a'^ are polynomials in and (3i are some combi- 
nations of Cj. The precise form of a^. Pi depends on the 
choice of the regulating function fA{l)- For a finite value of 
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the cut-off A the function Vcr. non-poi. in Eq. H2.11(l con- 
tains, in general, not only non-polynomial terms in q but 
also polynomial ones, which are however suppressed by 
inverse powers of the cut-off A. The terms in Eqs. H2.13(l . 
(|2.14|l proportional to A, A? and A? correspond to linear, 
quadratic and cubic divergences in Eqs. (|2.1() . (|2.2() and 
are absent in the dimensionally regularized expressions. 
Renormalization can be performed in very much the same 
way as before by absorbing the terms proportional to A^ 
A^ and Inyl by the counter terms. The only difference is 
that we now need the LO counter terms in order to get rid 
of the yl^-term in Eq. (|2.13|l at NLO and the yl^-term in 
Eq. H2.14|l at NNLO. In addition, NLO counter terms are 
required to renormalize the NNLO TPE potential. Notice 
that since Vcr, non-poi. depends on A, the renormalized 
expressions for the potential using dimensional and cut- 
off regularizations are only identical with each other for 
A —> oo. Does that mean that we should necessarily take 
A — > oo7 The answer is no. This is because in an EFT one 
is usually only able to calculate observables with a finite 
accuracy performing calculations up to a certain order in 
the low-momentum expansion. Taking A A^, the error 
from keeping A finite is beyond the theoretical accuracy. 
In other words, since 



Vcr , non — pol. . non— pol. 

+ 0{1/A), 



(2.15) 



the DR and CR expressions are identical up to higher 
order terms. At NNLO one should choose /a(0 in such a 

way that ycTno„-poL = ^DR?.on-pol. + 0{1/A^). 



2.2 Spectral function representation 

Before discussing implications of the choice of the regular- 
ization scheme on the convergence of the chiral expansion, 
we will show how cut-off regularization of the potential 
can be understood in terms of the spectral function rep- 
resentation. By their very nature, these spectral functions 
are the most natural objects to separate the long- and 
short -distance contributions to the NN potential in mo- 
mentum space. We will switch to the notation introduced 
in Ref.ini and express V^^'=' , yNNLO ^^g. 

+ Vriq) {<Ji • q) (0-2 • q) , 

yNNLO ^ ^^(^) ^ ^^(^) (^^ . (^^ . 

+ Wriq) (ri • T2) (ai • q) ((72 • q) . (2.16) 

The functions Vi{q) {Wi(q)) correspond to isoscalar (isovec- 
tor) parts of the potential and can in case of dimensional 
regularization (and also for A oo) be read off from 
Eqs. (Eiill, lEHl)- The subscripts C, S" and T stand for the 
central, spin-spin and tensor contributions, in order. 

The functions Vi{q) {Wi{q)) can be represented (mod- 
ulo terms polynomial in q'^)^ by a continuous superposi- 

^ In the present work we are only interested in the finite- 
range part of the two-pion exchange, which is given by the 
non-polynomial terms in momentum space. 



tion of Yukawa functions 
2 



V{q) 
W,iq) 



dfJ-fJ- 2 : 2 

J2M„ M + Q 



(2.17) 



where Pi{fJ-) and rii{fi) are the corresponding mass spec- 
tra (spectral functions). Note that subtracted dispersion 
integrals should be used in Eq. (|2.17|) if spectral func- 
tions do not decrease for large /i. This usually happens 
in the EFT calculations, where the spectral functions are 
obtained within the \ow-fj, expansion. Notice further that 
subtraction constants can be absorbed by the LECs corre- 
sponding to the short-range contact interactions and thus 
do not introduce any additional ambiguity. It is easy to 
see that pi{p) and ?7i(/i) can be obtained from Vi(q) and 
Wi{q) via: 



p.,{p) ^lm[V^{0+ -ip)\ , 
^,{p) = lm\w^{0+ - ifi) 



(2.18) 



These spectral functions contain the whole dynamics cor- 
responding to the exchanged tttt system. Once pt (/x) , iji (p,) 
functions are determined, the TPE potential can easily be 
obtained using Eq. 12.17|l . 

Let us now calculate the spectral function p^{p), which 
results from the integral in Eq. (|2.2|l regularized with a 
cut-off A. We will choose the regulating function /yi(0 ^ 
IaQ-) = d{A — I). Performing integration over angles, one 
obtains 



3^: 



di 



1287r2F4 7o g(/2 + ^2 + 4^/2) 
X (^8ci + C3 (/' - q^)^ 



(2.19) 



In (/ + g)2 + 4Af 2 - In (? - q)^ + 4A'l 



One then finds for the spectral function 

pi.{p)^lm[v^{0+ -tp) 
3.9^ 



64F4 
1 



(2M2(2ci - C3) + C3m') (2M2 - p^) 



X -9{p - d{^A^ + 4M2 - p) . (2.20) 

M 

The entire yl-dependence is contained in the second Heavi- 
side step-function, which represents the only difference to 
the DR expression. Thus, as one could expect for phys- 
ical reasons, cutting off the momentum I of one of the 
exchanged pions at ^ = yl leads to a cut-off in the TPE 
spectral functions, which in this specific case takes the 
value a/v12+4M^ ~ A. Clearly, the precise form of the 
resulting regulator in the spectral function representation 
depends on the choice of /yi(?). Similar relations between 
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the pion momentum and spectral function cut-ofFs can be 
obtained for other contributions to the NN potential as 
well. Using the regularized spectral function representa- 
tion for TPE opens therefore an easy and convenient way 
to obtain the CR expressions and will be adopted in what 
follows. To be specific, we define the CR spectral functions 
pf{fi), rif{ii) according to 



(2.21) 



where pi (/i) , rji {p) are the corresponding DR spectral func- 
tions, see section 12.31 The non-polynomial parts of the 
TPE potential at NLO and NNLO have then the same 
structure as in Eqs. (|2.4() . (|2.8(l . where the loop functions 
L{q) and A{q) should be replaced by L'^i.q) and A'^{q) 
defined as: 



L'^iq) ^e{A-2M^)—ln 



4M2(yl2-f g2) 



s = ^A^ - 4Af 2 , 



A%)-0iA^2AU)l.rct.nf^. (2.22) 

Several comments are in order. First of all we note that 
for 2M7r < A the CR and DR expressions Vcr and Vdr 
only differ from each other by higher order contact inter- 
actions (i.e. by short-range terms) if Vcr is expanded in 
powers of 1/A. One can therefore use this regularization 
prescription in calculations at any given order in the low- 
momentum expansion without getting into trouble with 
spurious long-range contributions suppressed by inverse 
powers of A which might arise for a different choice of the 
cut-off function. Further, one should keep in mind that 
our choice of regularization is by no means unique. Differ- 
ent choices lead to equivalent results for the potential (up 
to higher order terms) and may be used as well. Finally, 
we would like to point out that the spectral function rep- 
resentation (|2.17fl does not allow to properly reproduce 
terms, which are polynomial in q^ and non-analytic in 
M2. Such terms are not important for our present work 
since we do not consider variation in . If one is inter- 
ested in the M7r-dependence of the nuclear force, see for 
instance Ref. [32, a cut-off regularization should be per- 
formed at the level of divergent integrals in Eqs. H2.1|l and 
H2.2|l rather than in the spectral function representation. 



2.3 Coordinate space representation 

The coordinate space representations Vc,s,T{r) {Wc,s,T{r)) 
of the isoscalar (isovector) central, spin-spin and tensor 
parts of the potential Vc,s,T{q) {Wc,s,T{q)) are defined 
according to 

V{r) ^Vcir) + Wcir) (ti • T2) 

+ Ws{r) (ti • T2) ) (cti ■ (T2) 

+ {vT{r) + WT{r) (X1-T2)) 

X (3(71 • f cr2 • r - CTi • 0-2) • (2.23) 



The functions Vc',s,t('') can be obtained for any given 
r > from the corresponding spectral functions via 



Vc{r) = 
Vrir) = 

Vs{r) = 



1 



dppe ^''pc(m) , 



(2.24) 



2 A/, 



1 



dppe '"' {3 + 3pr + p'^r'^)pT{p) , 



2AU 



(2.25) 



2M, 



dppe {^p^pt{p) - 'ips{^^) 



(2.26) 



The coordinate space representation of the isovector parts 
is given by the above equations replacing Vc.s,T{'r) 

Wc,s,T{r) and pc,S,t{p) iic,s,t{p)- 

We will now consider the coordinate space representa- 
tion of the CR TPE potential. According to the definition 
(Einil, Eqs. (I23I), lEHl) and itO^ we obtain for the CR 
spectral functions at NLO 



(2.27) 



P 



3g\ ^p^ - 4A/2 



e{A - p) 



and at NNLO 



Pc(m) 



3gi 

64Fi 



(2M2(2ci - C3) 



C'iP 



X {2Ml - p^) - 0{p - 2AU) d{A - p) , 



4{p)^\4{p) 
p^ 



9l 



128FS 



Ci{AMl - p^) 



x-9{p-2M^)e{A- p) . 
P 



(2.28) 



We are now in the position to discuss implications 
of keeping the momentum space cut-off A in the above 
expressions for the spectral functions finite. Let us con- 
sider, for example, the isoscalar central part of the TPE. 
In Fig. 121 we show the (normalized) integrand I{p) in 
Eq. as a function of ^ for r = ^M'^ r = Af^\ 

r = 2M~^ and for yl 00. For the LECs ci 3 we use 
the values a = -0.81 GeY~\ C3 = -4.7 GeV-\ 
As expected, at large distances r > 2M^^ the integral 
in Eq. H2.24|l is dominated by low-p components in the 
spectrum, where the chiral expansion for the spectral func- 
tion is well behaved. However, at intermediate distances 
r 1/A/^ ~ 1.4 fm the dominant contribution to the in- 
tegral comes already from the region p 0.6 GeV, where 
only a very slow (if at all) convergence of the chiral ex- 
pansion for p{p) is expected. Certainly, at even shorter 
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r = 0.5M^"' 







(1 [GeV] 



Fig. 2. The (normalized) integrand /(/i) in 
Eq. 12.241 for different distances r. 



distances the resulting TPE potential is completely deter- 
mined by the region of large where the spectral function 
is not properly described in chiral EFT. It is then clear 
that setting the upper limit of the integral in Eq. (|2.24|l to 
infinity, which corresponds to the DR result, leads to inclu- 
sion of spurious short-range physics in the TPE potential. 
We stress again that this problem solely arises because at 
the order we are working there are no contact terms which 
normally would absorb these contributions. On the other 
hand, introducing a finite cut-off A in the spectral func- 
tion representation according to Eq. H2.21(l we explicitly 
exclude all short-range components (i.e. those ones with 
the range R < A^^) from the TPE potential, which are 
still present in the DR expression. The procedure is legit- 
imate and does not lead to any ambiguity if A is chosen 
to be of the order of (or larger than) A^^ ~ Mp, which 
provides a natural scale for the effective field theory. Var- 
ious choices for A lead to exactly the same result for low- 
energy observables provided one keeps terms in all orders 
in the EFT expansion. Choosing a specific value for A in 
the calculation at any finite order in the low-momentum 
expansion one implicitly makes a particular choice for the 
combination of the higher order contact terms. There are 
obviously no restrictions in choosing A if all LECs are of 
the natural size and the expansion parameter in EFT is 
small. Both dimensional and cut-off regularizations lead 
to similar results and the unphysical short-range compo- 
nents of TPE resulting from keeping A very large (or even 
oo) are compensated by corresponding counter terms to 
the order at which one is working. In fact, as we will show 
in the next section, the NN potential at NLO may serve 
as an example of such a situation. 

In certain cases, however, it appears to be advanta- 
geous to explicitly remove the spurious short-distance physics 
when calculating chiral loops. This happens, for instance, 
in SU(3) baryon chiral perturbation theory, where a much 
slower convergence of the chiral expansion is expected due 
to the relatively large mass of the strange quark. The lead- 
ing nonanalytic components from loop corrections calcu- 
lated using dimensional regularization in some cases seem 
to destroy the good agreement of the lowest order calcula- 



tion with data (like e.g. the Gell-Mann-Okubo relation for 
the baryon masses). Of course, this disagreement is cor- 
rected after inclusion of higher order contributions. How- 
ever, the chiral expansion seems not to behave well and no 
clear convergence can be observed to the orders yet calcu- 
lated (such a statement holds e.g. for the baryon masses 
but not for the baryon magnetic moments 34V5F). 
Reformulating chiral EFT using a finite cut-off regulariza- 
tion allows to remove the spurious short-distance physics 
and to improve the convergence |26ll27ll28j . if the cut-off 
procedure is implemented model-independently (as shown 
in Ref. 

As discussed in the introduction, slow convergence of 
the chiral series for the NN interaction is observed if one 
uses the values for the LECs 01,3^4 obtained in ttN scatter- 
ing. The numerically large values of 03^4 lead to a strong 
and attractive TPE at NNLO. Strong deviations from the 
NPSA results are observed in the D- and F-waves which 
are parameter-free at this order and are still sensitive to 
the TPE contribution. We will demonstrate in the next 
section that the problems with the convergence are not 
due to the large values of 03^4 which provide a proper long- 
range part of TPE |23II24| . but rather due to unphysical 
short-range components in the DR expressions for the po- 
tential, which can (and should) be avoided using the CR. 

To close this section we will give analytical expressions 
for TPE in the coordinate space. For the NNLO contribu- 
tions one finds: 



-2a; 



327r2Fi 



2\2 



2ci x^{l + x^) 



-Fc3(6-hl2a;-hl0a;2-f 



1287r2Fi 



+ C3(24 + y(24+12y + 4t/2 + y3) 



-Ax 



'■{2 + 2y + y^)+4x'') 



(2.29) 



9\ 



-2x 



487r2Fi 



-C4 (1 + x)(3 + 3x -I- x^) 



9\ 



7687r2Fi 



e-y ( 



48- 



24y2 + 



-4x^(8 + 5y + 2/2) 



(2.30) 



9l 



-2x 



487r2Fi 



-C4(l-l-x)(3-h3x-^2x2) 



9l 



3847r2F4 r" 
-4x^(2 + 2?/ 



— ^C4 24 + 24y + 12?;^ + Ay" 



(2.31) 



with the abbreviations x = M^r, y — Ar. For the NLO 
contributions we could not perform integrations in Eqs. I|2.24f) - 
(j2.26|l analytically. In the chiral limit {M^ = 0) the results 
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take, however, the foUowing simple form: 



M^=0 



23g\ - I0g\ - 1 
15367r3Kfr5 



- 6 



15 



9% 



M„=0 



+e-^(15 + 15y + 6y2 + y^) 
6 



9\ 



-e'y{6 + 6y + 3y^+y^) 



(2.32) 



(2.33) 



(2.34) 



Certainly, the large-distance asymptotical behavior of the 
potential is unaffected by the cut-off procedure provided 
that A ^ M^r. In Fig. Owe compare the isoscalar central 
part of TPE obtained using CR and DR. The strongest ef- 
fects of the cut-off are observed at intermediate (~ M~^) 
and smaller distances, where TPE becomes unphysically 
attractive if DR is used. In contrast, removing the large 
components in the mass spectrum of the TPE with the 
reasonably chosen cut-off A = 500 . . . 800 MeV greatly 
reduces this unphysical attraction and the resulting po- 
tential is of the same order in magnitude as the one ob- 
tained in phenomenological boson-exchange models. It re- 
mains to say that the CR expressions for the TPE are still 
not regular in the origin, although the short-distance be- 
havior is milder than in the case of DR and the leading 
singularities at r = are removed. For example, while 
Vc{r) cx 1/r^ V^{r) cx l/r^ for r ^ 0. The Schrodinger 
(or LS) equation has still to be regularized by introduction 
of an additional cut-off. 



3 Chiral TPE at NNLO and peripheral NN 
scattering 

In this section we will apply the OPE and TPE nuclear 
force at NNLO, calculated using the cut-off regularization 
as described in the previous section, to NN phase shifts 
with orbital angular momentum / > 2 and to mixing an- 
gles with J > 2. As already pointed out in the introduc- 
tion, no contact terms contribute at NNLO to the scat- 
tering amplitude in these channels. Consequently, such pe- 
ripheral phase shifts are entirely determined by the long- 
range part of the nuclear force and thus provide a sensitive 
test of the chiral TPE. The OPE potential at NLO (and 
NNLO) is given by 

V 9\ (, 4M2 \ (^i-q)(a2-q) 

I/qPE = --TE^ ^ "18 \Ti-T2 2 I ;i/f2 ' 

(3.1) 

where the LEC dig, is related to the Goldberger-Treiman 
discrepancy via 

^ = f^fl-^rfrs). (3.2) 



In what follows we use qa — 1-26, i^^ = 92.4 MeV, dis = 
-0.97 GeV-2, which leads to g^N ^ 13.2. For TPE at 
NLO and NNLO we use Eqs. ^T^, (|TSI) with the fu nction s 
L(g), A{q) replaced by i^(g), ^^(g) defined in Eq. ^TF^. 
For the c^'s we adopt the values [T51: ci = —0.81 GeV^"'^, 
C3 = -4.70 GeV-\ C4 = 3.40 GeV-i. 

The explicit expressions for the partial-wave decom- 
position of the potential can be found in Ref. |23. The 
partial-wave projected Lippmann-Schwinger (LS) equa- 
tion for the NN T~matrix reads: 



p2 _ (p, 



/l\2 



t;:^j,(p",p), (3.3) 



where the on-shell S- and T-matrices are related via 



Si I' [P, P) ^Sii'l 



87r2 



pmTi ],{p, p) . 



(3.4) 



Since the chiral potential grows at large momenta, the LS 
equation (|3.3|l has to be regularized, ^ which requires in- 
troduction of an additional cut-off (which can be done 
consistently with the CR discussed so far). It is com- 
monly believed (and observed in various boson-exchange 
and phenomenological potential models, see also Ref. [TF|) 
that because of the centrifugal barrier the NN interaction 
in the peripheral partial waves becomes weak enough to 
be treated perturbatively. This is also confirmed by the 
smallness of the corresponding phase shifts. To calculate 
phase shifts in high partial waves one may therefore use 
the Born approximation to the T-matrix:^ 



T[\,{p',p)^Vi%ip',p). 



(3.5) 



Such a procedure, which is analogous to the one of Refs. [S] 
155] . allows to avoid the introduction of an additional cut- 
off in the LS equation and will be adopted in the present 
work. One should, however, always keep in mind that this 
approximation breaks down if the phase shifts become 
large. 

Before presenting our predictions for the high partial 
waves, we would like to specify the differences between our 
analysis and the one of Kaiser et al. fSJ: 

— First of all, we strictly follow the lines of Weinberg's 
power counting and do not include relativistic correc- 
tions to the nuclear force at NNLO, which have been 
included in Ref. 9 . In contrast to this reference, we 
use non-relativistic kinematics when calculating phase 
shifts. 

— We also do not include the contribution from once it- 
erated OPE, which turns out to be numerically small 



^ This regularization of the LS equation should not be con- 
fused with the CR of the chiral loops discussed in previous 
sections. 

Notice, however, that the weakness of the NN interaction 
for high values of I is not related to the chiral expansion and 
does not follow from the power counting. 
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Fig. 3. The isoscalar central TPE potential at NNLO in r-space. The solid line shows the DR result 
corresponding to /l = cx) while the red band results from varying A between 500 and 800 MeV. The short- 
(long-) dashed line shows the phenomenological a {a + lo + p) contributions based on the isospin triplet 
configuration space version (OBEPR) of the Bonn potential EH- 



in most channels. From the point of view of the power 
counting, there is no reason to include once iterated 
OPE and not to include two-, three-, . . . times iter- 
ated OPE. 

— We use slightly different (and more modern) values for 
the LECs ci,3_4 and for gA- 

— The last and most important difference is that we use 
the finite cut-off A to regularize chiral loop integrals. 
This is in strong contrast with the analysis of Ref. 0| , 
where DR corresponding to yl = oo has been adopted. 
As it is clear from the above discussion, optimal val- 
ues for the cut-off are those close to the scale where 
the EFT description becomes inaccurate. Taking a too 
small A will remove the truly long-distance physics 
while too large values for the cut-off may affect the 
convergence of the EFT expansion due to inclusion of 
spurious short-distance physics. We will therefore vary 
A in the range A = 500 . . . 800 MeV which appears to 
be physically reasonable and matches well with both 
Mp and the cut-off used in the LS equation for NN 
scattering [TH] . 

Let us start with the D-waves which are shown in 
Fig. 21 The LO result represented by pure OPE already 
provides a good approximation to the phase shifts in the 
^Di and partial waves and to the mixing angle 62- 
It is too weak in the channel and does not describe 
properly the ■^£'3 phase shift. The letter appears to be 
quite small (|^| ~ 4.6° at Euh = 300 MeV) compared to 
the other D-wave phase shifts {\S\ ~ 9.7° . . .25.5°). The 
reason is that partial-wave projected OPE, taken on the 
energy shell, is strongly suppressed in this channel. Con- 
sequently, the ^D^ phase shift is quite sensitive to TPE 
but also to the iteration of the potential which we neglect 
in the present analysis. 



The NLO predictions obtained using dimensional reg- 
ularization are shown by the dashed curves. One observes 
a visible improvement for the ^Z?2 phase shift and for £2 
while the NLO corrections go in the wrong direction in 
the "^Di and "^D^ channels. The NNLO predictions calcu- 
lated with dimensional regularization are depicted by the 
solid lines. The good agreement with the data observed at 
LO and NLO is destroyed in all partial waves for energies 
Eiah > 50 MeV and the chiral expansion does not seem 
to converge. Note that at N'^LO one independent contact 
operator contributes to each D-wave so that the agree- 
ment with the data will presumably be restored, see also 
Ref. for a related discussion. Note that the results 
presented here are parameter-free and are very similar to 
the ones of Ref. (01 . No breakdown of the chiral expansion 
at NNLO is observed using the CR in the chiral loops and 
choosing A = 500 . . . 800 MeV. This proves explicitly that 
strong disagreement with the data resulting from chiral 
TPE at NNLO is due to unphysical short-distance com- 
ponents, which are kept in the DR expressions for the po- 
tential. The use of a momentum space cut-off keeps only 
the long-distance portion of the chiral loops and leads to 
a greatly improved behavior. Our predictions agree with 
the data in the ^I?2 and ■^Z?i channels and go into the right 
direction in the '^D^ partial wave and for €2- Note that the 
uncertainty due to the variation of A is quite significant 
at higher energies. The NLO bands are narrower than the 
NNLO ones as a consequence of the fact that the leading 
TPE contribution is numerically quite small. 

Our predictions for the F-wave phase shifts and for ea 
are presented in Fig. O Although the situation with DR 
subleading TPE is not as dramatic as for D-waves, a too 
strong attraction is clearly visible in the '^F2, ^F^ and '^^4 
partial waves. Removing the short-distance components of 
TPE with the cut-off regularization leads to a significant 
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Lab. Energy [MeV] 

Fig. 4. D-wave NN phase shifts and mixing angle 62 versus the nucleon laboratory energy. The dotted curve 
is the LO result (i.e. pure OPE), while the dashed (solid) curve refers to NLO (NNLO) results for OPE+TPE 
with the potential calculated using dimensional regularization. The light (dark) shaded band shows the NLO 
(NNLO) predictions with chiral TPE obtained using the cut-off regularization with A = 500 . . . 800 MeV. The 
filled circles depict the Nijmegen PSA results |37|. 



improvement in the and channels, while additional 
repulsion is still missing in the partial wave. 

For the G- and higher partial waves we observe very 
similar results for both DR and CR TPE and the phase 
shifts are essentially given by OPE. Notice that addi- 
tional attraction is generated from iterated OPE in the 
^Gs partial wave which is not included in the present 
work. Therefore, the strong disagreement with the data in 



this channel is presumably just an artefact of the Born- 
approximation. 

Finally, we would like to point out that the typical un- 
certainty of ^ 30% observed in our NNLO predictions at 
^-lab = 300 MeV is consistent with the power counting. 
Indeed, the N-^LO counter terms are expected to provide 
corrections to the S-matrix of the order Q'^/A^. Taking 
Q ~ 375 MeV, which corresponds to E^h = 300 MeV, 
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Fig. 5. F-wave NN phase shifts and mixing angle £3 versus the nucleon laboratory energy. For notation see 
Fig.El 



and identifying with the smallest value used for the 
cut-off A, i.e. A = 500 MeV, one estimates the N^LO ef- 
fects as ~ 32%. As already stated before, the uncertainty 
is larger in the cases, where phase shifts are numerically 
small (like, for instance, in the ^i?2, '^D^, partial 
waves). Although the following remark is quite obvious, 
we stress that a two-nucleon potential based on a system- 
atic EFT approach should not be fine-tuned to fulfill a 
X^/dof ~ 1 as it is done in more conventional approaches. 
This does, however, not mean that such a precision can 
not be reached. 



4 Summary and conclusions 

In this paper we have considered the two-nucleon poten- 
tial in chiral effective field theory, making use of a novel 
method of regularizing the pion loop integrals. For that, 
we have considered the spectral functions obtained from 
the NLO and NNLO TPE contributions and argued that 
only masses below the chiral symmetry breaking scale 
should contribute in the loop integrals. This can be easily 
implemented by applying a cut-off to the spectral func- 
tions. Varying the cut-off between 500 and 800 MeV (as 
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Lab. Energy [MeV] 



Fig. 6. G-wave NN phase shifts and mixing angle £4 versus the nucleon laboratory energy. For notation see 
Fig.El 



given by the mass of the heaviest Goldstone boson, the 
eta, and the mass of the lightest resonance, the rho), we 
find the following pertinent results: 

1) From the regularized spectral functions, we have con- 
structed the coordinate space representations of the 
various components of the NN potential. The isoscalar 
central TPE shown in Fig. 13 agrees with phenomeno- 
logical potentials. The strong unphysical attraction of 
the TPE is a short-range phenomenon which is sup- 
pressed in CR by choosing A < Mp, whereas in the 
Bonn potential np exchange cancels the strong attrac- 



tive TPE, which is one particular model for this kind 
of short-range physics. 
2) We have considered the peripheral partial waves (/ > 
2), because at NNLO, these are given entirely by OPE 
and TPE with no free parameters. We have calculated 
these phases in Born approximation which should be 
legitimate at least for the D- and higher waves, see 
Ref. |18| for more discussion. The uncertainty in most 
D- and F-waves at NNLO is still sizeable even with 
the finite cut-off. This has to be expected because of 
the large values of Ci 's. The results for D- and F-waves 
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are still not completely converged at NNLO, but the 
error of about 10 (1)° at E'lab = 300 MeV for the D- 
(F-)waves appears reasonable. There is no breakdown 
of the chiral expansion for D-waves beyond Tiab = 50 
MeV and for F-waves beyond Tiab = 150 MeV as found 
in Ref. [5] using dimensional regularization. 

3) It is no surprise that the NNLO TPE gives larger cor- 
rections than the NLO one because of the delta contri- 
butions subsumed in the LECs C3 and C4. One might 
therefore contemplate including the Z\(1232) explic- 
itly in the effective Lagrangian (see e.g. [7]), because 
in such a theory most (but not all!) of the NNLO ef- 
fects are shifted to NLO, provided that a systematic 
analysis of pion-nucleon scattering in such a scheme is 
available (for attempts see e.g. '39''40",'4I ). However, 
for obtaining a precise potential one still would have 
to go to NNLO, which is considerably more compli- 
cated than in the pion-nucleon EFT, as witnessed by 
the fact that no complete fourth order calculation with 
deltas in the single nucleon sector exists. 

4) We stress that dimensional regularization is by no means 
ruled out by such considerations. In general, for quickly 
converging expansions, it should be the method of choice. 
If, however, the convergence for some well understood 
physical reason is slow and (some) observables become 
sensitive to spurious short-distance physics kept in 
DR, it might be preferable to use CR, as done here. 
In our case the LECs C34 are large. Choosing DR, one 
generates a series of higher order contact interactions, 
which are oc C34 and therefore large. The D- and F- 
waves are certainly most sensitive to such contact in- 
teractions and are thus strongly affected. 

5) The cut-off used here is not a form-factor and is not 
related to the finite extension of the nucleon. Obvi- 
ously, the precise shape of the regulating function and 
the precise value of A are not important (as long as 
the value for A is well above the two-pion thresh- 
old, A > 2Mt^ and below the scale of chiral symmetry 
breaking, A < A^ ). 

In a subsequent publication, we will apply CR to the 
low partial waves in the non-perturbative regime, where 
we have to solve the regularized Lippmann-Schwinger equa- 
tion Eq. H3.3|) to generate the bound and scattering states. 
It can be demonstrated that there are no deeply bound 
states and low-energy observables are not affected by CR. 
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